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FIG. 1: Cu'' trimeric chains in Ca3Cu3(P04)4. The strongly 
coupled Cu'' trimer consists of a central square planar Cu(l) 
ion (black circle) and two pyramidal Cu(2) ions (gray circles). 
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The nuclear spin-lattice relaxation time Ti of '^^P nuclei in the title compound is measured for the 
first time and interpreted in terms of a modified spin-wave theory. We establish a novel scenario for 
one-dimensional ferrimagnetic spin dynamics — it is three-magnon processes enhanced by exchange 
scattering, rather than Raman processes, that make the major contribution to l/Ti. 

PACS numbers: 76.50.-l-g, 76.60.-k, 75.50.Gg 



I. INTRODUCTION 

Low-energy dynamics in one-dimensional (ID) quan- 
tum magnets is a long-standing problem and nuclear 
magnetic resonance (NMR) findings are often enlight- 
ening in this context. It was a fine collaboration 
that the NMR relaxation rates l/Ti and I/T2G (Rcf. 
0) of the spin-i antiferromagnet Sr2Cu03 illuminated 
multiplicative logarithmic corrections to the dynamic 
susceptibility^ in critical spin chains. Ti measurements 
on the spin-| antiferromagnet (CH3)4NMnCl3 were a pi- 
oneering study in an attempt to detect the long-time dif- 
fusive spin dynamics in one dimension^ The diffusive 
contribution to 1/Ti was observed in the spin-i quantum 
limit^'^ and for spin-gapped antiferromagnets^ as well. 

Antiferromagnets, whether critical or gapped, are thus 
vigorously studied, while little is known about ferri- 
magnetic, as well as ferromagnetic, dynamics. Re- 
cently numerous ID ferrimagnets have been synthesized 
in an effort to design molecule-based ferromagnets. The 
static properties of various heterospin chains were corre- 
spondingly calculatedfiSiSiifliiLiSiiiii but few effort has 
been devoted to exploring their dynamic features. Re- 
cent NMR observations of ID ferrimagnetsiSiiSiii have 
thus been motivated and have indeed stimulated the 
theoretical interest in them. The most pioneering Ti 
measurementsjA^ performed on the metal-radical hy- 
brid compound Mn(C5H202F6)2CioHi7N202, were in- 
terpreted in terms of solitonic excitations at low tem- 
peratures and based on the spin-diffusion model at high 
temperatures. However, this material has rather large 
exchange interactions and its magnetic susceptibility (x) 
times temperature (T) exhibits no ferrimagnetic mini- 
mum at a measurable temperature.^® The nonnegiigi- 
ble single-ion anisotropy and the spread magnetic mo- 
ment over the radical may also blur the intrinsic ferri- 
magnetic features. The family of manganeseporphyrin- 
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FIG. 2: (a) Dispersion relations of the spin-wave excitations, 
(b) Modified spin- wave calculations of the susceptibility com- 
pared with experimental findings sd, H — 0.1 T. 



based ferrimagnets, whose exchange interactions are 
somewhat smaller, was also observed through NMRjii 
but the complicated crystalline structure restricted the 
analysis to solving the magnetic configuration. 
NMR findingslSi on the bimetallic chain compound 
NiCu(C7H6N206)(H20)3-2H20 implied a spin-wave re- 
laxation scenario, but the averaging effect over the nu- 
merous protons masks the low-temperature dynamics. 

In such circumstances, we perform an NMR study on 
the homometallic chain compound Ca3Cu3(P04)4 (Fig. 

which is a novel ferrimagnet of topological origin.— 
This collaboration aims at verifying a newly developed 
modified spin- wave nuclear magnetic relaxation theory^ 
Two of the present authors have recently demonstrated 
that l/Ti in ID heterometallic ferrimagnets may sig- 
nificantly be enhanced by exchange interactions, in an 
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attempt to interpret the novel NMR observationaiSi 
for NiCu(C7H6N206)(H20)3-2H20, which are hardly 
understandable within the usual Raman relaxation 
scheme. However, all the findings^ were far from 
conclusive to a possibility of more than two magnons 
mediating the proton spin relaxation. Since only 
the nuclei were practically available as probes in 
NiCu(C7H6N206)(H20)3-2H20, the crystal water of nui- 
sance and the resultant ill-behaving spin-echo recovery 
curve with increasing field at low temperatures restricted 
the Ti analysis to rather high temperatures T ^ 60K. 
Ca3Cu3(P04)4 contains P atoms as efficient probes. The 
well-isolated Cu^^ chains with Ca columns in between, no 
single-ion anisotropy of the Cu^' spins, and the ferrimag- 
netic minimum of at a moderate temperature [see 
Fig. Olb)] furthermore guarantee this material to be a 
fine stage of ID ferrimagnetic dynamics. 



II. MULTIMAGNON-MEDIATED NUCLEAR 
SPIN RELAXATION 

We describe Ca3Cu3(P04)4 by the Hamiltonian, 



[see Fig. I3a)]. Ho gives two-body interactions and 
makes a crucial contribution to nuclear spin-lattice relax- 
ation. The 0{S~^) terms are neglected in the following. 
The dispersive branches ujk-.i and ujk-.s are reminiscent of 
the dual excitations in alternating-spin chains,^"' whereas 
the flat band ujk:2, describing intratrimer excitations, is 
peculiar to the present system. 

Our of modifying the conventional spin- wave 

theory is distinct from the original idea proposed by 
TakahashiS^ and Hirsch et al.m^ Their way of suppress- 
ing the divergent sublattice magnetizations consists of 
diagonalizing the Hamiltonian together with a Lagrange 
multiplier subject to zero staggered magnetization. The 
thus-obtained energy spectrum depends on temperature 
and fails, for instance, to reproduce the Schottky-peaked 
specific heat. Seeking after better thermodynamics, we 
diagonalize the bare Hamiltonian and then minimize the 
free energy with a Lagrange multiplier subject to zero 
staggered magnetization,^^ The thus-calculated x is in 
good agreement with observations [Fig. [2Jb)]. 

The hyperfine interaction between a '^^P nucleus 
and Cu^^ spins consists of isotropic Fermi contact and 
anisotropic dipolar coupling and is defined as 
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H. — [Jl(S'„:i • Sn:2 + Sn;2 ' Sn:3) 

n=l 

+ J2iSn+l:l ■ Sn;2 + S n;2 ' S'„_i:3) 

^gii^H{S:,,,+S';,,^ + Sl)\, (1) 

where each Cu(l) ion is antiferromagnetically coupled 
to four Cu(2) ions in an applied field {H) (see Fig. 

We set Ji/k-Q and J2/^b equal to 100 K and 
8K, respectivelyii^ Employing the Holstein-Primakoff 
transformationist we expand the Hamiltonian with re- 
spect to l/S* as 
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-2s^{Ji + J2)N + ni + no + o{s-^) 
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where Hi contains the 0(5"*) terms. Hi describes linear 
spin-wave excitations and is diagonalized in the momen- 
tum space as 

Hi = Ji + ./2) + ^ ( t^fc + ^ ujk:ial.,^ak:z I , (3) 

where a[..j creates a spin wave of ferromagnetic (i = 1, 2) 
or antiferromagnetic (i — 3) aspect, whose excitation en- 
ergy is given by 

S 

t^fcii = Wfc - — (Ji -I- J2) -l-.g/iB-ff, 

i^k:2 = {Ji + J2)S + gnBH, (4) 

^k:3 = ujk + —{Ji + J2) ~ gmH, 



with 



Wfc = I Y ( Ji + J2)2 + 32 Ji J2 sin^ ^, (5) 



Hm = g^iBh-fNl^ J2^'^^n:tSn:t + K-.t^n-.z)- (6) 

n i=l 

Ho and Hm arc both much smaller than Hi and may be 
regarded as perturbations to the linear spin- wave system. 
When we calculate up to second order in V = T^o + Hm, 
the probability of a nuclear spin being scattered from the 
state oi — m to that of = m -I- 1 is given by 
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Ei — E„ 



m-Ef), 



(7) 

where i and / denote the initial and final states of the 
unperturbed electronic-nuclear spin system. Then we ob- 
tain Ti = {I—m){I+m+l)/2W. Considering the signifi- 
cant difference between the nuclear and electronic energy 
scales at moderate fields and assuming the Fourier com- 
ponents of the coupling tensors to have little momentum 
dependence as Y.n e"^("+^/2-i)^A _ ^ ^ ^ ^j^^, ^la.- 
man and three-magnon relaxation rates read 
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FIG. 3: (Color online) Various nuclear spin-lattice relaxation processes. Spin waves which are emitted in the first-order 
mechanism (solid arrows) flip a nuclear spin (x) via the hyperfine interaction, whereas four-magnon exchange correlations 
(dotted arrows) thermally scatter a first-order virtual spin wave with q = — ^4, where ferromagnetic spin waves are drawn by 
black (ujk-.i) and colored {ujk:2) straight arrows, while antiferromagnetic ones (LOk-.a) by wavy arrows, (a) First-order Raman 
processes, (b) First- and second-order three-magnon processes, which are inseparable in nonlinear equations. 
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where n^-.i = {<xl..^ak:i) is the thermal distribution func- 
tion of modified spin waves, while k2 ^ and fcg*'* are deter- 



mined through LJki : 
ujki-.i ~ Tt-^-n = 0, 



0, UJi^(l),^+UJk2:l 



— TiwN 

"3 

and Wj,(3).3 - Wfej:! - Wfc^:3 -I- huj^ = 0. Wii{ki,k2) and 

Wy7(fci, fe, fcs) are diagrammatically represented in Fig. 
13 Within the first-order mechanism containing a nuclear 
spin in direct contact with spin waves via the hyperfine 
interaction, any multi-magnon relaxation rate is much 
smaller than the Raman one. However, some of multi- 
magnon processes make a significant contribution to 1 /Ti 
through the second-order mechanism, where a nuclear 
spin flips with the help of virtual spin waves which are 
then scattered thermally via the four-magnon exchange 
interaction. We consider the leading second-order relax- 
ation, that is, exchange-scattering-induced three-magnon 
processes, as well as the first-order relaxation. Second- 
order single-magnon and Raman relaxation processes, 
containing three and two virtual magnons, respectively, 
are much more accidental due to the momentum con- 
servation and much less contributive due to the magnon 
series damping. As for four-magnon processes, the first- 
order relaxation is nonexistent to begin with, whereas the 
second-order one originates in the six-magnon exchange 
interaction and contains two virtual magnons. Thus and 
thus, all other higher-order processes have no significant 
effect on the relaxation scenario. 



III. EXPERIMENTAL TEST 

Now we present ^"'^P NMR spectra for powder samples 
of Ca3Cu3(P04)4 in Fig. 01 The crystallographically 
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FIG. 4: NMR spectra measured at 51.711 MHz. They 
can be decomposed into typical powder patterns [the insets 
(a) and (b)], originating in the P(l) (dotted lines) and P(2) 
(broken lines) sites, and the NMR shift for the latter is plotted 
as a function of the susceptibility at H — 3T [the inset (c)] . 



inequivalent P sites, labeled P(l) and P(2), give two dis- 
tinct lines, each of which has a characteristic shape of 
the anisotropic powder pattern, as is revealed in the in- 
sets (a) and (b). Both lines broaden and shift to lower 
field with decreasing temperature, implying that contact 
as well as dipolar terms exist in the hyperfine field on 
the P nuclei. The NMR shifts are thus anisotropic and 
their parallel (K^^) and perpendicular {K±) components 
for the P(2) sites are plotted as functions of the suscepti- 
bility at = 3 T in the inset (c). If we roughly estimate 
the hyperfine coupling between each P(2) nucleus and 
the nearest-neighbor Cu^^ ion through the relations A± = 
NAtJ.BdK±/dx and ^|| = A^AMBdi^n/dx, where A^a is the 
Avogadro number, we find that A± = 0.66kOe//iB and 
All = 1.58kOe//iB. In order to reveal the low-frequency 
spin dynamics, we measure Ti for the P(2) sites through 
the saturation-recovery method. 

Temperature dependences of 1/Ti are shown in Fig. 
|S1 In each Cu^^ trimer unit, the two spins Sn-.i and 
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FIG. 5: (Color online) Experimental (symbols) and theoreti- 
cal (lines) findings of l/Ti as functions of temperature. I/t'^^ 
and l/rf are also plotted in (b). 
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FIG. 6: (a) Semilog plots of l/Ti as functions of l/T. Dot- 
ted Lines are guides to the probable exponential behavior 



1/Ti oc e"^off/^ at each field, (b) The thus-extracted Af^ 
as a function of field. The two slopes, Z\[,g = gusH/kB and 



l(3) 



2g^BH/kB, are sho-wn for reference. 



Sn:2 are almost equidistant from the nearby P(2) atom, 
"whereas the other one Sn;3 is much more distant from 
thati^ Therefore, the isotropic coupUng constants are 



taken as 



^ A- 



and =0. Since dipolar 



interactions are also sensitive to the location of corre- 
lating moments, the anisotropic coupling constants may 
be taken similarly as A\ — A\ = and A\ = 0. 
Then we set A~ and A^ equal to 0.85kOe//iB and 
1.20kOe//iB, respectively, which are both consistent well 
with the experimental findings A^ — 0.66kOe//iB and 
v4j| — 1.58kOe//iB. Considering that recent electron- 
spin-resonance measurements of this compound^^ have 
yield temperature-dependent and anisotropic g values 
(.9|| > 5± — 2), the value may be closer to the the- 

oretical parametrization. We are further convinced of the 
coupling constants employed finding the nearest Cu(l)- 
P(2) distance to be about (l/A^)!/^ = 2.5 A, which is 
in excellent agreement with the crystalline structure.SS, 
The thus-calculated l/Ti = 1/rf ^ + 1/Tf ^ reproduces 
the observations pretty well. The cxchange-scattering- 
enhanced three-magnon relaxation rate generally grows 
into a major contribution to 1/Ti with increasing tem- 

perature and decreasing field. While both 1/T{ and 
1/T-l exhibit an exponential behavior at low tempera- 
tures, their activation energies, referred to as fcB^cff ^'^'^ 
/cB^cff ' respectively, look different. At moderately low 



FIG. 7; Experimental (symbols) and theoretical (lines) find- 
ings of l/Ti as functions of field at 95 K (a) and 170 K (b), 
where and l/T^^ are also plotted. The observations 

are replotted as functions of 1/ \/H and fitted to an expression 
1/Ti =P/x/H + Q (c). 



temperatures and weak fields, fioj^ <C k^T <C J2, Eq. 
© reads 



1 



Ji + J; 



t[^) 2TThSJ_ 

(10) 

where Kq is the modified Bessel function of the sec- 
ond kind and behaves as Ko{hLLj-[^ /2kBT) ~ 0.80908 — 

(2) 

\n{huJN / ksT) . Thus we learn that A^^^ ~ gfi^H /k^. 
Equation © is much less analyzable, but Fig. I^b) 
claims that < 2i|,^\ Figure brings the next lead- 



ing exponential behavior l/Ti cx e '^Itt/^ to light, be- 
cause Eq. (Cni) is vahd for T < Ja/fcB 8K. A^^^ looks 
like 2g/iB-ff/fcB or more, rather than g^L^H /k^. The two- 
magnon-mediated nuclear spin relaxation is mainly given 
by H^ii(fci,fc2), where a ferromagnetic spin wave of en- 
ergy oJki-i is created and that of energy ujk2:i — '^ki-.i 
is destructed, while the three-magnon-mediated one by 
Wiii{ki,k2,k'i), where two net spin waves of energy 
Wfej:! and Wfcj:! are created and that of energy tok^-i ~ 
ujk-i-.i +'^fe3:i is destructed. Therefore, Raman processes 
are activated by the energy ujk=a:i — g^^BH, which is 
consistent with Eq. (|10|) . whereas three-magnon ones 
roughly by the energy 2a;fc=o:i — '^gi^BH, which may look 
somewhat larger due to the complicated field dependence 
lying in the second-order mechanism. 

Field dependences of l/Ti at higher temperatures 
more impress on us the significance of three-magnon 
processes. Figures |7{a) and |7{b) show that the ac- 
celerated relaxation with decreasing field can never be 
explained by the Raman scheme but should be at- 
tributed to exchange-scattering-enhanced three-magnon 
processes. The spin-diffusion model may be mentioned 
in this context. Diffusion-dominated ID spin dynam- 
ics gives l/Ti of the form P/^/H + Qy^ where the first 
and second terms come from transverse and longitudi- 
nal spin fluctuations, respectively, and are both posi- 
tive. Figure [7{c) shows that the present observations 
fitted to the diffusive law result in negative Q. We do 
not exclude a possibility of diffusive dynamics appearing 
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in ID ferrimagnets as well, but a distinct field depen- 
dence of the second-order relaxation mechanism masks 
such a moderate field effect in the present case. It was 
in AgVP2S6 rather than in the most familiar Haldane- 
gap antiferromagnet Ni(C2H8N2)2N02C104, whose exci- 
tation spectrum drastically varies with increasing field, 
that integral-spin diffusive correlations were observed^ 
Exchange-scattering-induced three-magnon processes are 
sensitive to an excitation gap and their contribution to 
1 /Ti is strongly suppressed, for example, by slight mag- 
netic anisotropy. There are indeed some indications of 
spin diffusion— in the ferrimagnctic chain compound 
Mn(C5H202F6)2CioHi7N202 with nonncgligiblc single- 
ion anisotropyiiS 

IV. CONCLUDING REMARKS 

We have performed NMR measurements on the topo- 
logical fcrrimagnet Ca3Cu3(P04)4 and have confirmed 
a novel scenario for ID spin dynamics^ multimagnon- 
mediated nuclear spin relaxation, by showing the 
parametrization to be crystallographically convincing, re- 
vealing that the relaxation is activated by twice the gap 
rather than the gap itself at low temperatures, while it is 
remarkably accelerated with decreasing field, and point- 
ing out the irrelevance of the spin-diffusion model. In- 
deed pioneering Ti observations on the layered ferromag- 



net CrCls (Ref. ^ and the coupled-chain antiferromag- 
net CsMnCls -2112 (Ref. i32t} claimed to have detected 
three-magnon processes, but they were both, in some 
sense, classical findings under the 3D long-range order^ 
Our findings are literally ID quantum spin relaxation 
beyond the Raman mechanism, which were obtained 
through an elaborately modified spin-wave theorytS 

The bond-alternating homometallic chain compound 
Cu(C5ClH4N)2(N3)2 (Ref. H is another anisotropy-free 
ferrimagnet of topological origin^^ and therefore NMR 
measurements on it are highly encouraged. Since ID fer- 
romagnets, which can be regarded as low-energy sectors 
of ID ferrimagnets, may also play this fascinating sce- 
nario, more understanding will come with further exper- 
iments, for instance, on the spin-i ferromagnetic chain 
compound (CH3)4NCuCl3ifi& The quasi-one-dimensional 
mixed-spin ferromagnet MnNi(N02)4(C2ll8N2)2 (Ref. 
Isri l is also highly interesting in this context, whose low- 
energy spectrum consists of two dispersive ferromagnetic 
excitation branches 
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